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PURPOSE: To efficiently calculate an optimal design 
value by properly setting variable difference amount In 
the design of lens, eto., by Ihe method of attenuation 
least square. 

CONSTITUTION: A CPU 1 calculates the eigen value 
distribution information on Ihe product matrix (A^TAv) 
between the Jacoblan matrix Av havir>g the variable 
amount of itie charBCteristic (evaluation function pv (Xv)) 
of the system for the variable difference amount gXj of 
the system being design object as an element and 
transpose of matrix AvT and sets the variable difference 
amount 5Xj by using the eigen value distribution 
Infonnation. Instead of the eigen value of the product 
matrix (AvTAv), the square value of the singular value 
of the Jacobian matrix A may be used. 
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(54) [Title of the Invention] METHOD AND DEVICE FOR 
SUPPORTING DESIGN OF LENS OR THE LIKE 

(57) [Abstract] 

[Object] To efficiently calculate an optimal design value 
by adequately setting a variable difference amount in the 
design of a lens or the like by an attenuation least square 
method. 

[Solving Means] A CPU 1 calculates eigen value 
distribution information on the product matrix (Av'^Av) of the 
Jacobian matrix Av with the variable amount of the 
characteristic (an evaluation function Fv(Xv) ) of the system 
for the variable difference amount 8x3 of the system being a 
design object as elements thereof and the transposed matrix 
Av*^ and sets the variable difference eimount 6xj by using the 
eigen value distribution information. The square value of 
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the singular value of the Jacobian matrix A may be used in 
place of the eigen value of the product matrix (Av'^Av) • 
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[ Claims ] 

[Claim 1] A design supporting device to support the design 
of a lens or the like by an attenuation least square method 
for locally changing and optimizing a plurality of 
components so that a system of a design object comprising 
the plurality of components has a desired characteristic 
comprising : 

an eigen value distribution information operating means 
to operate eigen value distribution information on the 
product matrix (Av'^Av) constituted of the product of the 
Jacobian matrix (Av) with the variable amount of the 
characteristic of the system to the difference amount of the 
variables corresponding to the plurality of components as 
the elements thereof and the transposed matrix (Av'^) of the 
Jacobian matrix; and 

a difference amount setting means to set the difference 
amount of the variables by using the eigen value 
distribution information . 

[Claim 2] The design supporting device according to Claim 1, 
wherein the difference amount of the variables is set so 
that the initial value of the damping factor in the 
attenuation least square method is present between the 
minimum value and the maximum value of the eigen value 
distribution of the product matrix. 

[Claim 3] The design supporting device according to Claim 1, 
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wherein the difference amount of the variables is set so 
that the median of the eigen value distribution of the 
product matrix is substantially matched with the initial 
value of the damping factor in the attenuation least square 
method . 

[Claim 4] The design supporting device according to Claim 1, 
wherein the difference amount of the variables is set so 
that the mean of the eigen value distribution of the product 
matrix is substantially matched with the initial value of 
the damping factor in the attenuation least square method. 
[Claim 5] The design supporting device according to Claim 1, 
wherein the difference amount of the variables is set by 
using eigen value distribution information on the product 
matrix and predetermined weighting information set 
corresponding to the eigen value distribution. 
[Claim 6] A design supporting device to support the design 
of a lens or the like by an attenuation least square method 
for locally changing and optimizing a plurality of 
components so that a system of a design object comprising 
the plurality of components has a desired characteristic 
comprising : 

a square value distribution information operating means 
to operate square value distribution information on the 
eigen value of the Jacobian matrix (Av) with the variable 
amount of the characteristic of the system to the difference 
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amount of the variables corresponding to the plurality of 
components as the elements thereof; and 

a difference amount setting means to set the difference 
amount of the variables by using the square value 
distribution information . 

[Claim 7] The design supporting device according to Claim 6, 
wherein the difference amount of the variables is set so 
that the initial value of the damping factor in the 
attenuation least square method is present between the 
minimum value and the maximum value of the square value 
distribution of the eigen values of the Jacobian matrix. 
[Claim 8] The design supporting device according to Claim 6, 
wherein the difference amount of the variables is set so 
that the median of the square value distribution of the 
eigen value of the Jacobian matrix is substantially matched 
with the initial value of the damping factor in the 
attenuation least square method. 

[Claim 9] The design supporting device according to Claim 6, 
wherein the difference amount of the variables is set so 
that the mean of the square value distribution of the eigen 
value of the Jacobian matrix is substantially matched with 
the initial value of the damping factor in the attenuation 
least square method, 

[Claim 10] The design supporting device according to Claim 
6, wherein the difference amount of the variables is set by 
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using square value distribution information on the eigen 
value of the Jacobian matrix and predetermined weighting 
information set corresponding to the square value 
distribution . 

[Claim 11] A design supporting method to support the design 
of a lens or the like by an attenuation least square method 
for locally changing and optimizing a plurality of 
components so that a system of a design object comprising 
the plurality of components has a desired characteristic 
comprising : 

a step of operating eigen value distribution 
information on the product matrix (Av'^Av) constituted of the 
product of the Jacobian matrix (Av) with the variable amount 
of the characteristic of the system to the difference amount 
of the variables corresponding to the plurality of 
components as the elements thereof and the transposed matrix 
(Av*^) of the Jacobian matrix; and 

a step of setting the difference amount of the 
variables by using the eigen value distribution information. 
[Claim 12] The design supporting method according to Claim 
11, wherein the difference amount of the variables is set so 
that the initial value of the damping factor in the 
attenuation least square method is present between the 
minimum value and the maximum value of the eigen value 
distribution of the product matrix. 
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[Claim 13] The design supporting method according to Claim 
11, wherein the difference eunount of the variables is set so 
that the median of the eigen value distribution of the 
product matrix is substantially matched with the initial 
value of the damping factor in the attenuation least square 
method, 

[Claim 14] The design supporting method according to Claim 
11, wherein the difference amount of the variables is set so 
that the mean of the eigen value distribution of the product 
matrix is substantially matched with the initial value of 
the damping factor in the attenuation least square method. 
[Claim 15] The design supporting method according to Claim 
11, wherein the difference amount of the variables is set by 
using eigen value distribution information on the product 
matrix and predetermined weighting information set 
corresponding to the eigen value distribution. 
[Claim 16] A design supporting method to support the design 
of a lens or the like by an attenuation least square method 
for locally changing and optimizing a plurality of 
components so that a system of a design object comprising 
the plurality of components has a desired characteristic 
comprising : 

a step of operating square value distribution 
information on the eigen value of the Jacobian matrix (Av) 
with the variable amount of the characteristic of the system 
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to the difference amount of the variables corresponding to 
the plurality of components as the elements thereof; and 

a step of setting the difference amount of the 
variables by using the square value distribution information. 
[Claim 17] The design supporting method according to Claim 
16, wherein the difference amount of the variables is set so 
that the initial value of the damping factor in the 
attenuation least square method is present between the 
minimum value and the maximum value of the square value 
distribution of the eigen values of the Jacobian matrix. 
[Claim 18] The design supporting method according to Claim 
16, wherein the difference amount of the variables is set so 
that the median of the square value distribution of the 
eigen value of the Jacobian matrix is substantially matched 
with the initial value of the damping factor in the 
attenuation least square method. 

[Claim 19] The design supporting method according to Claim 
16, wherein the difference amount of the variables is set so 
that the mean of the square value distribution of the eigen 
value of the Jacobian matrix is substantially matched with 
the initial value of the damping factor in the attenuation 
least square method. 

[Claim 20] The design supporting method according to Claim 
16, wherein the difference amount of the variables is set by 
using square value distribution information on the eigen 
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value of the Jacobian matrix and predetermined weighting 
information set corresponding to the square value 
distribution . 

Detailed Description of the Invention 
[0001] 

[Technical Field of the Invention] The present invention 
relates to a design supporting device and a design 
supporting method using an attenuation least square method 
out of optimization methods to locally changing a plurality 
of components so that an objective system comprising the 
plurality of components have desired characteristics, and 
for example, it relates to a design supporting device and a 
design supporting method to support a lens design, an 
optimum pattern design of an LSI, a building design and the 
like . 
[0002] 

[Description of the Related Art] Components to constitute a 

system as a design object are expressed by variable vectors 

consisting of n pieces of elements as follows. 

[0003] 

[ Formula 1 ] 

Xv = [Xi X,]^ (1) 

And the amount to evaluate the characteristics of the 
object is expressed by an evaluation function vector 
consisting of m pieces of elements as follows. 
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[0004] 
[Formula 2] 

Fv(Xv) = [Fi(Xv) F^CXJ]"^ (2) 

[0005] 

[ Formula 3 ] 

F,,(X^) = w,,{f^(X^)-f^.,,,} (3) 

Where, suffix v is added to indicate Xv and Fv are 
vectors (or matrix). Further, fk(Xv) , fk and tar are the 
values of the evaluation function to indicate respective 
characteristic values, and a target value thereof, and wk is 
a weighting coefficient. A constraint value may be included 
in these evaluation functions in addition to the 
characteristic value intended for a desired value. In the 
attenuation least square method, a solution is generally 
obtained so that the amount of deviation (error) from a 
target value of each evaluation function value is minimized, 
and an objective function expressed by the following formula 
is used for a single evaluation scale. 
[0006] 
[ Formula 4 ] 

(l)(Xv) = F/(Xv)Fv(Xv) (4) 

The objective function expressed by the formula (4) is 
referred to as a merit function in the field of the lens 
design, in particular. By applying the least square method 
to the formula (4), the necessary information of Xv for a 
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local minimum value is given as follows . 

[0007] 

[Formula 5] 

(A/A^)AXv = -A/Fvo (5) 

Where, Av is a Jacobian (Jacobi) matrix of Fv(Xv) , and 
defined by Formula (6). 
[0008] 


[ Formula 6 ] 



d Fi 

d Fi 


d XI 

d X2 


d Fz 

a F2 


d XI 

d X2 


5Fm 

• 

aFm 


d XI 

d X2 


(6) 


d Fi 

d Xn 

a Fz 

d Xn 
i 

d Xn 


In Formula (6), Fi means an i-th element Fi(Xv) of the 
evaluation function vector Fv(Xv) (Formula (2)). AXv is a 
solution vector to indicate fluctuation of the variable 
vector to the next step, and expressed as follows. 
[0009] 
[Formula 7] 

AX^ = [AX, AXJ"^ (7) 

Generally, in the attenuation least square method, an 
optimum solution is obtained by solving Formula (8) by using 
a parcimeter (a damping factor) p to correct non- linear 
components of the secondary or higher differential in place 
of directly solving Formula (5). 
[0010] 
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[Formula 8] 


vo 


(8) 


Where, Iv is a unit matrix. 
[0011] In the field of the lens design, a term of the 
attenuation least square method (DLS method: a Damped Least 
Squares method) is extensively used- On the other hand, in 
the field of the non-linear planning method or the like, a 
term of the Levenberg-Marquardt method is extensively used. 
[0012] On the other hand, in the actual numerical 
calculation, each element a^j of Av is calculated not by 
using a strict differential coefficient value (Formula (9)) 
expressed by Formula (6), but calculated by using Formula 
(10) for the difference value. 
[0013] 
[ Formula 9 ] 


[0014] 

[Problems to be Solved by the Invention] However, the value 
of the optimum variable difference amount 8xj is diversified 
according to the characteristic of the object system, 
raising a problem in that, if an inadequate value is given, 
errors occur in the numerical computation, many iterative, 
calculations are required, and the optimum solution cannot 
be efficiently obtained. 


a^j = aF^/aXj 
= 8Fi/8xj 


(9) 


(10) 
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[0015] The present invention is achieved in light of this 
point, and an object of the present invention is to provide 
a design supporting device and a design supporting method 
capable of obtaining the optimum solution in a constantly 
efficient manner by adequately setting the variable 
difference amount. 
[0016] 

[Means for Solving the Problems] In order to achieve the 
above -described object of the present invention, a design 
supporting device to support the design of a lens or the 
like by an attenuation least square method for locally 
changing and optimizing a plurality of components so that a 
system of a design object comprising the plurality of 
components has a desired characteristic comprises an eigen 
value distribution information operating means to operate 
eigen value distribution information on the product matrix 
(Av'^Av) constituted of the product of the Jacobian matrix 
(Av) with the variable amount of the characteristic of a 
system to the difference amount of the variables 
corresponding to the plurality of components as the elements 
thereof and the transposed matrix (Av*^) of the Jacobian 
matrix, and a difference amount setting means to set the 
difference amount of the variables by using the eigen value 
distribution information . 

[0017] In order to achieve the same object of the present 
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invention, in a design supporting method to support the 
design of a lens or the like by an attenuation least square 
method for locally changing and optimizing a plurality of 
components so that a system of a design object comprising 
the plurality of components has a desired characteristic, 
eigen value distribution information is operated on the 
product matrix (Av'^Av) constituted of the product of the 
Jacobian matrix (Av) with the variable amount of the 
characteristic of the system to the difference amount of the 
variables corresponding to the plurality of components as 
the elements thereof and the transposed matrix (Av*^) of the 
Jacobian matrix, and the difference amount of the variables 
is set by using the eigen value distribution information. 
[0018] Further, preferably in the design supporting device 
or the design supporting method, the difference amount of 
the variables is set so that the initial value of the 
damping factor in the attenuation least square method is 
present between the minimum value and the maximum value of 
the eigen value distribution of the product matrix, the 
difference amount of the variables is set so that the median 
of the eigen value distribution of the product matrix is 
substantially matched with the initial value of the damping 
factor in the attenuation least square method, the 
difference amount of the variables is set so that the mean 
of the eigen value distribution of the product matrix is 
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substantially matched with the initial value of the damping 
factor in the attenuation least square method, and the 
difference amount of the variables is set by using eigen 
value distribution information on the product matrix and 
predetermined weighting information set corresponding to the 
eigen value distribution. 

[0019] In addition, in the design supporting device or the 
design supporting method, the "square value distribution of 
the singular values of the Jacobian matrix (Av)" may be used 

in place of the "eigen value distribution of the 
product matrix (Av'^Av)". 
[0020] 

[Operation] According to the device of Claim 1 or according 
to the method of Claim 11, eigen value distribution 
information on the product matrix (Av'^Av) of the Jacobian 
matrix (Av) with the variable amount of the characteristic 
of the system to the difference amount of the variables 
corresponding to a plurality of components as the elements 
thereof and the transposed matrix (Av'*') of the Jacobian 
matrix, and the difference amount of the variables is set by 
using the eigen value distribution information. 
[0021] According to the device of Claim 6 or according to 
the method of Claim 16, square value distribution 
information on the singular values of the Jacobian matrix 
(Av) with the variable amount of the characteristic of the 
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system to the difference amount of the variables 
corresponding to a plurality of components as the elements 
thereof is operated, and the difference amount of the 
variables is set by using the square value distribution 
information . 
[0022] 

[Embodiments] Embodiments of the present invention will be 
described with reference to the drawings . 

[0023] Fig. 1 shows a configuration of a design supporting 
device according to an embodiment of the present invention. 
This device comprises a CPU (a Central Processing Unit) 1 to 
perform various kinds of operations to optimize a variable 
corresponding to a component of a system for the design 
object, an input device 2 to input data and operation 
commands to the CPU 1 by a user, a storage device 3 to store 
programs executed by the CPU 1 and data or the like under 
operations, and an output device 4 to display the result of 
operations on a display device or to print the result of 
operations by a printer. These components 1 to 4 are 
connected to each other by a bus 5 . 

[0024] Next, the outline of the processing executed by the 
CPU 1 will be described. 

[002 5] Firstly, the product matrix (Av'^Av) of the above- 
described Jacobian matrix Av and a transposed matrix Av^ 
thereof is subjected to the eigen value break-down as 
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follows . 


[0026] 


[Formula 10] 


(11) 


Where, Vv denotes an orthogonal matrix, and Sv denotes 
an diagonal matrix of n, n-type with n-eigen values of the 
product matrix (Av^Av) as elements as shown in Formula (12). 
[0027] 

[Formula 11] 


The distribution of n- eigen values si, ... sn which 
are thus obtained is generally changed as the difference 
cimount 8xj of the above -de scribed variable is changed. In 
this case, the following method is employed in the present 
embodiment to directly control n-eigen value distribution by 
the difference eimount 6x j . 

[0028] Firstly, attention is paid to a linear approximation 
in a vicinity of XvO of the above -described evaluation 
function vector Fv(Xv) as expressed below. 
[0029] 

[Formula 12] 


Si 


0 


Sv = 


(12) 


0 


Fv(Xv) = Fv(Xvo) + A^AX^ 


(13) 


Then, the i-th element of the evaluation function 
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vector Fv{Xv) will be expressed as follows (Formula (14)) 
[0030] 

[Formula 13] 

3 F. 5 F. 

Fi(X,} = Fi(Xp) + (-g^)AXl-f • • ■+(-y^)AX„ 


A X. 


AX. 


5 X. 


In Formula (14), 8jFi denotes the variable amount by the 
i-th variable of the i-th evaluation function. From the 
relationship of Formula (14), it is proved that completely 
similar optimum operation can be performed even by using the 
matrix Av' and the vector AXv' such as in Formulae (15) and 
(16) in place of the Jacobian matrix Av (Formula (6)) and 
the solution vector AXv (Formula (7)). 
[0031] 

[Formula 14] 


Av'- 


6 iFl 62F1 (JnFl 

6 1F2 5 2F2 • • • 6 nFz 


S iFm 5 2F11 


(5 nFir 


(15) 


(16) 


By this method, the value of each element of the 
Jacobian matrix can be directly changed by the difference 
amount 8xj of the variable, and n-eigen values of si to sn 
can be freely controlled. 
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[0032] As described above, it is determined whether or not 
the set variable difference amount 6xj is an adequate value 
from the relationship between distribution of n-eigen values 
si to sn by the difference amount 6xj of the variable which 
is once given and the initial value of the damping factor p, 
and the optimum difference amount is set by changing the 
variable difference amount 8xj so as to obtain a more 
adequate value by using the information, and the optimal 
design value is efficiently obtained, 

[0033] Next, the present invention will be described more 
specifically by an example of the lens design. 
[0034] Fig. 2 shows a thin-walled lens model as the design 
object according to the first embodiment of the present 
invention. In Fig. 2, the lens has a thickness so as to 
easily show the embodiment. However, the lens thickness and 
the spacing in the calculation are set to be zero. The 
radius of curvature r3 is used so that the focal distance of 
the entire system is correctly 1 to this lens system, and 
remaining two radii of curvature rl and r2 are defined as 
variables. The tertiary spherical difference coefficient 
and the coma aberration coefficient are employed as two 
evaluation functions. The initial shape is expressed by the 
coordinates of (l/rl, l/r2) as follows. 

1/rl = -1.8, and l/r2 = 0.5 (17) 

The difference amounts 6(l/rl) and 6(l/r2) of these two 
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variables (1/rl) and (l/r2) are set in a range specified as 

follows . 

[0035] 

[ Formula 1 5 ] 

8{l/ri) = 1 X lO'^-l X 10\ 8(l/rJ = 1 x lO'^'-l x 10° (18) 
In addition, the variable difference amount to be set 
may be selected from four cases for each variable, i.e., 
sixteen in total at every 10^ in the range of Formula (18). 
Generally, in the lens design, the range of the variable 
difference amount to be roughly set is empirically 
understood according to the kind of the variables. 
[0036] In this condition, the Jacobian matrix Av for two 
variables and two evaluation functions in each variable 
difference amount set value is obtained. Two eigen values 
si and s2 of the product matrix (Av'^Av) are calculated as 
follows (entirely expressed by Formula (19)). 
(Margin below) 
[0037] 

[Formula 16] 

6(l/rJ = 1 X 10'', ad/rj = 1 x 10"' 

Si = 8,6xl0"'\ S2 = 2.2 X lO'MFig. 3(a)) 

8(l/rJ = 1 X 10-', 8(l/r2) = 1 x 10"' 

Si = 8,0 X 10"^°, S2 = 2.3 X 10"' (Fig. 3(b) ) 

6(l/ri) = 1 X 10"', 8(l/r2) = 1 x 10"' 

Si = 9.0 X 10'^°, S2 = 2.1 X 10"^ (Fig. 3(c) ) 
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6(l/rJ = 

= 1 

X 

10" 

*. 8(l/r,) 

= 

1 

X 

10" 



= 5 

.5 

X 

10-', s, = 

6 . 

8 

X 

10° (Fig. 

3(d) ) 

6(l/rJ = 

= 1 

X 

10- 

\ 6{l/rJ 

= 

1 

X 

10"* 



= 8 

.6 

X 

10-". = 

2 

. 2 

X 

10-* (Fig. 

4(a) 

6(l/rJ = 

= 1 

X 

10- 

8(l/r,) 

= 

1 

X 

lo- 



= 8 

.6 

X 

10-', s, = 

2 . 

2 

X 

10"" (Fig. 

4(b)) 

6(l/r,) = 

= 1 

X 

10" 

\ 8{l/r,) 

= 

1 

X 

10"^ 



= 8 

.2 

X 

10-'. = 

2 . 

3 

X 

10"" (Fig. 

4(c) ) 

6(l/rJ = 

= 1 

X 

10" 


= 

1 

X 

10° 


Si 

= 5 

.5 

X 

lo-^ = 

6. 

8 

X 

10° (Fig. 

4(d) ) 

6(l/r,) = 

= 1 

X 

10" 

\ 8(l/rJ 

= 

1 

X 

10"* 


Si 

= 8 

.6 

X 

10-". = 

2 

. 2 

X 

10° (Fig. 

5(a)) 

8(l/rJ = 

= 1 

X 

10" 

\ 6(l/rJ 

= 

1 

X 

10"- 


Si 

= 8 

.6 

X 

10-'. s, = 

2. 

2 

X 

10° (Fig. 

5(b) ) 

6(l/r,) = 

= 1 

X 

10" 

\ 6(l/rJ 

= 

1 

X 

10-^ 



= 2 

.2 

X 

10". = 

8.- 

7 : 

X : 

L0-' (Fig. 

5(c) ) 

6(l/rJ = 

= 1 

X 

10" 

\ 6{l/r,) 

= 

1 

X 

10" 



= 4 

.3 

X 

lo-^ = 

8 . 

5 

X 

10° (Fig. 

5(d) ) 

6(l/rJ = 

= 1 

X 

10" 

, 8(l/r,) 

= 

1 

X 

10-* 


Si 

= 5 

.8 

X 

10-", S3 = 

5 

. 2 

X 

10"' (Fig. 

. 6(a) 
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10" 



1 

X 
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Si 

= 5 

.8 
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10-', S3 = 

5. 
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10*' (Fig. 

6(b) ) 

6(l/r,) = 

= 1 
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10" 

, 8(l/rJ 
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10*' (Fig. 

6(c)) 

8(l/rJ = 

= 1 

X 

10" 

. 8(l/r3) 


1 

X 

10° 



- 22 - 


Si = 1.5 X 10\ S2 = 5-2 X 10"' (Fig, 6(d) ) 

(19) 

As described above, the eigen value distribution of 
each product matrix (Av'^Av) is indicated in Fig. 3(a) to Fig. 
6(d). These figures indicate the eigen value distribution 
operated by the CPU 1 and printed by the output device 4 . 
In each figure, the numerals are exponentially indicated on 
the horizontal axis, and the position of each eigen value 
(si, s2) is indicated by *. Further, the position of the 
initial value pO of the damping factor p is indicated by X 
on the horizontal axis, and the search range of the damping 
factor p therearound (hereinafter, referred to as "p search 
range") is indicated by [##... ##]. In the present 
embodiment , pO is set as follows . 
pO = 1.0 (20) 

In addition, the p search range is set as follows. 
[0038] 

[Formula 17] 

pO X 10'^ - pO X 0.5 X 10® (21) 

In the present invention, an adequate value of the 
difference amount of the variable is set from eigen value 
distribution information on the product matrix (Av'^Av) , and 
it is empirically understood in the field of the lens design 
to set one including the eigen values si and s2 in the p 
search range as many as possible. In this case, as clearly 
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shown in Figs. 3 to 6 , one including no eigen values si and 
s2 in the p search range is a case shown below (Formula 
(22) ) . 
[0039] 

[Formula 18] 

8(l/rJ = 1 X 10 *. 8(l/r2) = 1 x lo ' (Fig. 3(a)) 
8(l/rJ = 1 X 10 ', bil/r^) = 1 X lO" (Fig. 3(b)) 

(22) 

Further, one including one eigen value (si or s2) is 
the case shown below (Formula (23)). 
[0040] 

[Formula 19] 
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= 1 

X 
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(Fig. 

6(b) ) 


(23) 

Further, one including both the eigen values si and s2 
is the case shown below (Formula (24)), 
[0041] 

[Formula 20] 
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(Fig. 

6(d) ) 


(24) 

By setting any one of the variable difference amounts 
shown in Formula (24), iterative calculation to search the 
value of the damping factor p can be efficiently performed 
on the average. In order to set an adequate value of the 
variable difference amount more clearly, one in which the 
initial value pO of the damping factor p is smaller than the 
maximum value and larger than the minimum value of the eigen 
value distribution of the product matrix (Av'^Av) . 
[0042] 


8(l/rJ = 

1 

X 

10*. 

8(l/rJ 

= 1 

X 

10° 

(Fig. 

3(d) ) 

8(l/rJ = 

1 

X 

10*. 

8(l/rJ 

= 1 

X 

10° 

(Fig. 

4(d) ) 

6(l/r,) = 

1 

X 

10\ 

8(l/rJ 

= 1 

X 

10* 

(Fig. 

5(a)) 

6(l/rJ = 

1 

X 

io-\ 

8(l/rJ 

= 1 

X 

10" 

(Fig. 

5(b)) 

8(l/rJ = 

1 

X 

10 ^ 

8(l/rJ 

= 1 
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(25) 

In the present embodiment, this requirement is 
satisfied. 
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[0043] In addition, an adequate value of the variable 
difference amount is quantitatively set, one may be set, in 
which the initial value of the damping factor corresponding 
to a plurality of variables rl and r2 is substantially 
matched with the median or the mean of the eigen value 
distribution of the product matrix (Av'^Av) . Substantial 
matching is different by the object, and it is empirically 
understood that matching by 10^ to 10^ is sufficient for the 
case of the lens design. In other words, when the median or 
the mean of the eigen value distribution of the product 
matrix (Av'^Av) is expressed by sm, the difference amount 
6{l/rl) and 8{l/r2) of the variable set at that time if pO 
satisfies the relationship shown below. 
[0044] 

[Formula 22] 

s, X lO-'-s, X 10"' ^ Po ^ s, X lO'-s, X 10' (26) 
In the present embodiment, there are two variables, and 
the median and the mean of the eigen value distribution of 
the product matrix (Av'^Av) are matched with each other. When 
checking the mean of the eigen value distribution 
corresponding to the variable difference amount, it is 
proved that the relationship of Formula (26) is satisfied if 
being set like Formula (25). 

[0045] To set the variable difference amount more 
automatically, weighting as shown in Fig. 7 may be performed 
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so that the evaluation value becomes higher as more eigen 
values si and s2 are included in the p search range. In Fig. 
7 , the value on the horizontal axis indicates each eigen 
value, and the value on the vertical axis indicates the 
weight on each eigen value. If the eigen value is included 
in the p search range, 1 is set, otherwise 0 is set. The 
variable difference amount may be set so that the evaluation 
value is much higher thereby. 

[0046] Further, if the variable difference amounts 8(l/rl) 
and 8(l/r2) are set based on the relationship between the 
initial value pO of the damping factor and the eigen value 
distribution, continuous weighting may be performed so that 
the value becomes maximum at the initial value pO of the 
damping factor shown in Fig. 8, and the value is gradually 
decreased as being farther therefrom. In addition, 
weighting information (the weighting function) shown in Figs. 
7 and 8 is stored in the storage device 3. 

[0047] The result of the optimum operation by setting the 
above -described variable difference amount will be described 
below. 

[0048] Firstly, the number of the eigen values in the p 
search range, the mean of the merit function obtained by the 
result of calculation, and the mean of the iteration number 
required in the calculation will be indicated below. 
[0049] The merit function value is normalized by the value 


at the starting point. 


[0050] 


Number of eigen 
value 


Mean of merit 
function 


Mean of iteration 
number 


0 


0 . 713 


13 


1 


0.036 


10 


2 


0.034 


6.7 


As clearly shown from this result, the more the number 
of the eigen values included in the p search range is , the 
smaller the obtained merit function value is, and the 
smaller the iteration number is required. Therefore, the 
design values of the optimum rl and r2 can be obtained by 
setting the variable difference amount so that two eigen 
values are included in the p search range . 

[0051] Next, the result will be shown below for the case in 
which the initial value pO of the damping factor is present 
between the minimum value and the maximum value of the eigen 
value distribution and the case in which the initial value 
pO of the damping factor is not present therebetween. 
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[0052] 


Mean of merit Mean of iteration 
function number 


Present: 0.073 6.6 

Not present: 0,143 10.3 

As clearly shown from this result, by setting the 
variable difference amount so that the initial value pO of 
the damping factor is present between the minimum value and 
the maximum value of the eigen value distribution, the 
obtained merit function value becomes small, and the 
iteration number also becomes small, and the optimal design 
can be efficiently performed. 

[0053] In addition, the result will be shown below for the 
case in which the difference between the initial value pO of 
the damping factor and the mean of the eigen value 
distribution is expressed by the exponent difference, and 
the difference is 10°, 10\ 10*, 10^ and 10®, respectively. 
[0054] 

Exponent Mean of merit Mean of iteration 

difference function number 

0 0.066 6.6 

3 0.005 9.0 

4 0.033 13.0 

7 0.469 13.0 

8 0.958 13.0 
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As clearly shown from this result, the merit function 
value is minimized when the exponent difference is about 3, 
and the iteration number becomes less as the exponent 
difference is smaller. Therefore, the optimal design can be 
efficiently performed by setting the variable difference 
amount so that the exponent difference is within a range of 
2 and 3. 

[0055] Fig. 8 shows a configuration of a lens system as a 
design object according to the second embodiment of the 
present invention . 

[0056] In the present embodiment, twenty two variables 
consisting of lens components are used, and the beeim 
aberration is mainly used for the evaluation function to 
indicate the characteristic of the lens. The break- down of 
these twenty- two variables includes eleven surface radii r 
of curvature, and eleven kinds of spacing g, and the 
variable difference amounts are set to an optimum value in 
terms of the inverse 1/r of the radius of curvature in a 
range as follows. 
[0057] 

[Formula 23] 

6(l/r) = 1 X 10"*"1 X 10° (27) 

The variable difference amounts are set to an optimum 
value in terms of the spacing g in a range as follows . 
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[0058] 

[Formula 24] 

8(g) = 1 X 10'^"1 X 10° (28) 

When each variable difference amount is considered at 
the interval of 10^, all possible combinations will be given 
as follows . 
[0059] 

[Formula 25] 

7^' X 4'' « 8.3 X 10^^ cases 

Even if all these combinations are checked in 10"* 
second per combination, it would take an immense amount of 
time of about 263 years. This is one of the most difficult 
problems (a so-called NP-completion type problem) to obtain 
a strict solution in the combination-optimization problem. 
A genetic algorithm (Genetic Algorithm) has been extensively 
known as a method to efficiently solve the approximate 
solution of this kind of problem. The genetic algorithm was 
devised by J. H. Holland in 1960s, and this is a method for 
efficiently solving an optimization problem by simulating 
the essence of the biological evolution as it is on the 
computer . 

[0060] More specifically, this method is executed by the 
procedure as shown in Fig. 12. Firstly, in Step SI, the 
problem is modeled, and the gene type is set. For example, 
a numerical string (10100110) is used for the gene type. 


Next, a large number of individuals of different genes 
(numerical strings in this case) are formed (for example, 
10110001, 01101010, 11101100 ..,) (StepS2). And, the 
adaptability of the formed individuals is evaluated ( Step 
S3). In other words, each individual (series) is evaluated 
by using the predetermined evaluation function, and the 
individuals (series) of low evaluation will be selected in 
Step S4 • 

[0061] In Step S5, the not-selected individuals are 
multiplied to compensate the number of selected individuals , 
and the crossing-over (Step S6) and the mutation (Step S7) 
are performed. If the adaptability of the obtained 
individuals is insufficient. Steps S3 to S7 will be repeated. 
Here, the crossing-over is performed by switching specified 
parts by selecting specified gene pairs, and the mutation is 
performed by changing some parts of the gene at some 
probability. 

[0062] When the above -described processing is applied to 
the setting of the optimum variable difference cimount , a 
point is how the gene type is set and how the adaptability 
is evaluated. 

[0063] In the present embodiment, firstly, the series in 
which an arbitrary pair of difference eunount set values are 
expressed by the exponent are employed. In other words, the 
first to eleventh number of the series from the left to the 
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right are set as the exponent value of the difference amount 
on each radius r of curvature, the twelfth to twenty- second 
number of the series are set to be the exponent value of the 
difference amount on each spacing g, all the difference 
amounts of 1/r are set to be 1 x 10'^, and all the difference 
amounts of g are set to be 1 x 10'^, the gene row to indicate 
this state is expressed as follows. 

-6-6-6-6-6-6-6-6-6-6-6-3-3-3-3-3-3-3-3-3-3-3 (29) 
[0064] On the other hand, information on the eigen value 
distribution of the product matrix (Av'^Av) is used for the 
evaluation of the adaptability. More specifically, the 
result of higher evaluation will be obtained when the 
initial value pO of the damping factor is smaller than the 
maximum of the eigen values of the product matrix (Av'^Av) , 
and larger than the minimum thereof. 

[0065] When the adaptability is evaluated more strictly, 
the initial value pO is compared with the median or the mean 
of the eigen value distribution, and higher evaluation will 
be obtained when the initial value is substantially matched 
with the median or the mean of the eigen value distribution. 
[0066] In the present embodiment, weighting of the 
distribution as shown in Fig. 10 is performed, and the 
evaluation value becomes higher when the eigen value is not 
outside the p search range but within the range, and higher 
evaluation is obtained when the eigen value is close to the 
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initial value pO of the damping factor though it is within 
the p search range . 

[0067] Description will be given by a specifically 
numerical example. Firstly, in the variable difference 
amount , the radius r of curvature and the spacing g of the 
plane are set to be common as follows . 
[0068] 

[Formula 26] 

6(l/r) = 1 X 10"^ 

6(g) = 1 X 10' (30) 

In this condition, the eigen value distribution is 
examined as shown in Fig. 11(a). This state forms a 
starting point, and the weighting distribution in Fig. 10 is 
multiplied by the eigen value distribution in Fig. 11(a) to 
obtain the evaluation value, and processing is performed to 
set an adequate variable difference cimount by the genetic 
algorithm. The variable difference amount set value 
obtained as the result is as shown below for the radius r of 
curvature and the spacing g, respectively. 
[0069] 

[Formula 27] 

8(l/r) = 1 X 10-^-1 x 10"^ 

8(g) = 1 X lO'-l X 10' (31) 

The eigen value distribution corresponding to the 

variable difference amount set value (31) is shown in Fig. 
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11(b). As clearly shown from the comparison between Fig. 
11(a) and Fig. 11(b), it is shown that the eigen value 
distribution is gathered around the initial value pO of the 
damping factor pO (= 1.0) by the above-described processing. 
[0070] Comparison of the normalized merit function value 
and the iteration number required for the calculation before 
and after the above-described processing is as follows. 
[0071] 

Merit function Iteration number 
value of calculation 

Before processing 5.0 x 10'^ 3 

After processing 1.3 x 10"^ 3 

As clearly shown from this result, according to the 
present embodiment, the convergence efficiency of the 
attenuation least square method can be enhanced. 
[0072] Next, the outline of the processing executed by the 
CPU 1 in the third embodiment of the present invention will 
be described. 

[0073] Firstly, the Jacobian matrix Av is subjected to the 

eigen value break-down. 

[0074] 

[Formula 28] 

Av = UvPvVv'' (32) 

Where, UvVv denote orthonormal matrixes of m, m type 


and n, n type, respectively, and when m > n, Pv is expressed 

as follows. 

[0075] 

[Formula 29] 

PI c 


Pv = 


0 


pn 


(33) 


When m ^ n, Pv is expressed as follows, 
[0076] 

[Formula 30] 

pi 0 


Pv= 


Pm 


(34) 


Pv is a singular value matrix with n- or m- singular 
values of the Jacobian matrix Av as the elements thereof. 
The product Pv'^Pv of the singular value matrix is expressed 
as follows . 
[0077] 

[Formula 31] 


Pv'Pv = 


(35) 


PiT 


The distribution of the square values pl^, ... pn^ of 
the thus -obtained n- singular values is changed by the 
difference amount 8xj of the above -described variables. In 


this condition, in order to directly control the 
distribution of the square values pl^, ... pn^ of the thus- 
obtained n- singular values by the difference amount 8x j , a 
method (Formulae (13) to (16)) similar to that of the first 
embodiment is employed also in the present embodiment . 
Therefore, the square values pl^, ... pn^ of the thus- 
obtained n- singular values can be freely controlled by the 
difference amount 8x j . 

[0078] As described above, it is determined whether or not 
the set variable difference amount is adequate from the 
relationship of the square value distribution of the 
distribution of n- singular values of the once- given 
variables by the difference amount 8xj and the initial value 
of the damping factor p, and optimum difference amount is 
set by successively changing the variable difference amount 
6x3 so as to obtain a further adequate value by using the 
information so that the optimal design value is efficiently 
obtained. 

[0079] The above -described method is applied to the same 
lens system identical to that of the first embodiment, and 
the square values pl^ and k2^ of the singular values of the 
Jacobian matrix Av corresponding to Formula (19) are 
calculated as follows (the entire values are expressed by 
Formula ( 36 ) ) . 
[0080] 
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[Formula 32] 


8(l/rJ 


1 

X 

10'\ 

8(l/r,) 

= 1 X 10"* 





8 

.6x10 

P.^ = 

2.2x10"* (Fig. 3(a) ) 

8(l/rJ 

= 

1x10*, 8(l/rJ = 

1x10-* 





8 

. 0x10 

' \ P2 = 

2.3x10-' (Fig. 3(b) ) 

8(l/r,) 

= 

1 

X 

10'\ 

8(l/r,) 

= 1 X 10"'' 



Px^ 


9 

.0 X 

10 '\ p,' 

= 2 . 1 X 10"' (Fig 

. 3(c)) 

8(l/rJ 

= 

1 

X 

10 

SCl/r^) 

= 1 X 10° 





5 

. 5 X 

10-', Pa" 

= 6.8 X 10° (Fig. 

3(cl) ) 

6(l/rJ 

= 

1 

X 

10-*, 

8(l/rJ 

= 1 X 10* 



P.^ 


8 

. 6 X 

10-", p," 

= 2.2 X 10* (Fig 

. 4(a)) 

6(l/rJ 

= 

1 

X 

10 \ 

8(l/rJ 

= 1 X 10* 



P.^ 


8 

. 6 X 

10 P3^ 

= 2.2 X 10* (Fig. 

4(b) ) 

6{l/rJ 

= 

1 

X 

10-*, 

8(l/rJ 

= 1 X 10-" 



Px^ 


8 

.2 X 

10 *, p,' 

= 2.3 X 10"^ (Fig. 

4(c) ) 

6(l/rJ 

= 

1 

X 

io-\ 

6(l/rJ 

= 1 X 10° 



Px^ 


5 

. 5 X 

10-*, p," 

= 6.8 X 10° (Fig. 

4(d) ) 

8(l/r,) 

= 

1 

X 


6(l/r,) 

= 1 X 10* 



Px^ 


8 

.6 X 

10-", p.^ 

= 2.2 X 10° (Fig. 

5(a) ) 

6(l/rJ 


1 

X 

io-\ 

8(l/r,) 

= 1 X 10* 



Px^ 


8 

. 6 X 

10-', p." 

= 2.2 X 10° (Fig. 

5(b) ) 

8(l/rJ 


1 

X 

lo-^ 

8(l/r,) 

= 1 X 10'" 



Px^ 


2 

. 2 X 

10°, p." = 

= 8.7 X 10-^ (Fig. 

5(c) ) 

8{l/rJ 


1 

X 

10-% 

8(l/r,) 

= 1 X 10° 



Px^ 


4 

. 3 X 

io-\ p,^ 

= 8.5 X 10° (Fig. 

5(d) ) 
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6(l/rJ = 1 X 1Q\ 8(l/r2) = 1 x 10*' 

Pi^ = 5.8 X 10"'\ p^^ = 5.2 X lO^MFig. 6(a)) 
8(l/rJ = 1 X 10°, 8(l/r2) = 1 x lO^" 

Pi' = 5.8 X 10"', p/ = 5.2 X 10*MFig. 6(b)) 
6(l/rJ = 1 X 10°, 6(l/r2) = 1 x lO'" 

Pi' = 5.9 X 10"\ P2'' = 5.2 X lO'MFig. 6(c)) 
5(l/rJ = 1 X 10°, 8(l/r2) = 1 x 10° 

Pi^ = 1.5 X 10°, P2^ = 5.2 X lO'MFig. 6(d)) 

(36) 

As clearly shown above, the distribution of the square 
values of the singular values of the Jacobian matrix Av is 
totally same as the distribution of the eigen values of the 
product matrix (Av^Av) . In other words, the square values 
pl^ and p2^ of the singular values of the Jacobian matrix Av 
are equal to the eigen values of the product matrix (Av'^Av) . 
[0081] Therefore, in the first embodiment, replacement of 
the singular values of the eigen values si and s2 of the 
product matrix (Av'^Av) by the square values pl^ and p2^ 
corresponds to the present embodiment, and similar effect to 
the first embodiment can also be demonstrated by the present 
embodiment . 

[0082] Further, similar replacement can also be performed 

even in the second embodiment . 

[0083] 

[Advantages] As described above in detail, according to the 


design supporting device of Claim 1, or according to the 
design supporting method of Claim 11, eigen value 
distribution information on the product matrix (Av^Av) 
constituted of the product of the Jacobian matrix (Av) with 
the variable amount of the characteristic of the system to 
the difference amount of the variables corresponding to a 
plurality of components as the elements thereof and the 
transposed matrix (Av'^) of the Jacobian matrix is operated, 
and the difference amount of the variables is set by using 
the eigen value distribution information- Therefore, the 
optimum value by the descent method is efficiently operated, 
and the optimal design value can be efficiently obtained. 
[0084] Further, according to the design supporting device 
of Claim 6 or according to the design supporting method of 
Claim 16, square value distribution information on the 
singular values of the Jacobian matrix (Av) with the 
variable amount of the characteristic of the system to the 
difference amount of the variables corresponding to a 
plurality of components as the elements thereof, and the 
difference amount of the variables is set by using the 
square value distribution information, and similar effect 
can be demonstrated thereby. 
[Brief Description of the Drawings] 

[Fig. 1] Fig. 1 shows a configuration of a design 
supporting device according to an embodiment of the present 
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invention. 

[Fig. 2] Fig. 2 shows a configuration of a first lens 
system as a design object. 

[Fig. 3] Fig. 3 shows distribution of the eigen value (a 
square value of a singular value of a Jacobian matrix) of a 
product matrix of the Jacobian matrix and a transposed 
matrix thereof . 

[Fig. 4] Fig. 4 shows distribution of the eigen value (a 
square value of a singular value of a Jacobian matrix) of a 
product matrix of the Jacobian matrix and a transposed 
matrix thereof . 

[Fig. 5] Fig. 5 shows distribution of the eigen value (a 
square value of a singular value of a Jacobian matrix) of a 
product matrix of the Jacobian matrix and a transposed 
matrix thereof . 

[Fig. 6] Fig. 6 shows distribution of the eigen value (a 
square value of a singular value of a Jacobian matrix) a 
product matrix of the Jacobian matrix and a transposed 
matrix thereof . 

[Fig. 7] Fig. 7 shows an example of weighting information 
(a weighting function) . 

[Fig. 8] Fig. 8 shows an example of weighting information 
(a weighting function) . 

[Fig. 9] Fig. 9 shows a configuration of a second lens 
system as a design object. 


[Fig, 10] Fig. 10 shows an example of weighting information 
(a weighting function). 

[Fig. 11] Fig. 11 shows distribution of the eigen value (a 
square value of a singular value of a Jacobian matrix) a 
product matrix of the Jacobian matrix and a transposed 
matrix thereof . 

[Fig. 12] Fig. 12 shows a flowchart to indicate the 
procedure of a genetic algorithm. 
[Reference Numerals] 
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3 NUMERICAL VALUE 
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FIG. 12 
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